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C ! Abstract 

■ A search for Z boson decays to pairs of 5-quark jets has been performed in the fuU 
^SJ . dataset collected with the CDF detector at the Tevatron pp collider. After the selec- 
tion of a pure sample of bb events by means of the identification of secondary vertices 

^ ' from 6-quark decays, we have used two kinematic variables to further discriminate the 

electroweak bb production from QCD processes, and sought evidence for the Z decay 

■ in the dijet invariant mass distribution. An absolute background prediction allows 
. the extraction of an excess of events inconsistent with the background predictions by 

3.23cr but in good agreement with the amount and characteristics of the expected 
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0^ I signal. We then fit the mass distribution with an unbinncd likelihood technique, and 

^ ' obtain a. Z ^ bb signal amounting to 91 ± 30 ± 19 events. 

CD ■ 
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• 1—1 . 

X : 1 Introduction 



Z decays to 6-quark pairs are not exactly a unknown piece of Physics. Since 1992 the LEP 
experiments have detected several millions of them, and more have come from the polarized 
beams of the SLC The process is thus very well understood; the Z is one of the best known 
particles, and there can be no surprise in the thereabouts. At a proton-antiproton collider 
this particular process has never been seen before, though. The UA2 collaboration published 
in 1987 an analysis of jet data where they could spot the combined signal of W and Z decays 
to dijet pairsQ, but the decay of the Z to 6 quarks was not separated from the other hadronic 
decays. 

The Z decay to b quarks is the closest observable process to the expected decay of the 
Higgs boson, and — in view of the chances of a Higgs discovery in run 2 at the Tevatron — 
the understanding of the process, the knowledge of the expected mass resolution for a recon- 
structed decay, and the confidence with the kinematical tools that may help extracting it, 

^Phys. Lett. 186B (1987), 452. The search was later updated with a larger data sample which produced 
a signal of 5, 367 ± 958 events: Z. Phys. C49 (1991), 17. 
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are topics worth careful investigation. Moreover, a clean Z peak in a dijet mass distribution 
is useful for jet calibration, a relevant issue for the top quark mass measurement. 

In the following we present a search for the Z —>■ bb process in L = llOpb^^ of pp 
collisions collected by CDF during the years 1992-1995. We will show that the signal can be 
extracted by means of a very stringent selection that allows a reduction of the background 
by four orders of magnitude, while retaining in the final dataset a handful of Z events. 

This paper is organized as follows: in Section 2 we describe the CDF detector, the 
datasets used in the analysis, and the selection that leads to our final sample; in Section 3 
the counting experiment is described, and the results are discussed; Section 4 deals with the 
unbinned likelihood fits to the mass distribution. In Section 5 we present our conclusions. 

2 Data Selection 

The CDF detector has been described in detail elsewhere^. We only mention briefiy here 
those detector components most relevant to this analysis. The silicon vertex detector (SVX) 
and the central tracking chamber (CTC) are immersed in a 1.4 T axial field and provide 
the tracking and momentum measurement of charged particles. The SVX consists of four 
layers of silicon microstrip detectors and provides spatial measurement in the r — planeQ 
with a resolution of 15 fim^. The CTC is a cylindrical drift chamber containing 84 layers 
grouped into 9 alternating superlayers of axial and stereo wires. The calorimeters, divided 
into electromagnetic and hadronic components, cover the pseudorapidity range \t]\ < 4.2 and 
provide the jet energy measurement; they are segmented into projective towers subtending 
about 0.1 X 15° in ?7 — space. Central muon candidates are identified in two sets of muon 
chambers, located outside the calorimeters, as stubs extrapolating to charged tracks inside 
the solenoid. The data were collected by a three-level trigger system; the first two levels 
are provided by hardware modules, the third consists in software algorithms optimized for 
speed. 

At the Tevatron the Z production cross section has been measured both in the e+e~ and 
in the fi~^fi~ final states, and found to be az x BR{Z — > e+e~) = 0.235 ± 0.003 (stat.) ± 
0.005 (syst.) ± 0.020 (lum.) nb and az x BR{Z /x+zx") = 0.202 ± 0.016 {stat.) ± 0.020 
(syst.) ± 0.017 (lum.) ntQ. By scaling the e~^e~ figures up to the b quark branching fraction^ 
one expects that about 110,000 Z decays to 6-quark pairs have taken place at CDF during 

2 F.Abe et al, Nucl. Instr. Meth. Phys. Res. A 271 (1988), 387, and references therein. 

•^In CDF the positive z axis lies along the proton direction, r is the radius from the z axis, 9 is the polar 
angle, and is the azimuthal angle. 

D.Amidei et al, Nucl. Instr. Meth. Phys. Res. A 350 (1994), 73; S.Tkaczyk, Nucl. Instr. Meth. 
Phys. Res. A 368 (1995), 179, and references therein. 

^ See for instance P.Quintas, Proc. XI Symposium on Hadron Collider Physics, World Scientific, Singa- 
pore 1996. 

^ We can use for that purpose the world average branching fractions T ^^j^i/T z^hadrons ~ 0.2212±0.0019 
and Tz-.hadronslTz^e+e- = 20.77 ± 0.08 (PDG 1996). 
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run 1, the cross section times branching ratio for the process resulting in 1.080±0.029(stat.© 
syst.) ± 0.092(/um.) nanobarns. 

The natural trigger for a Z — > 66 signal would be a low energy jet trigger. However, the 
rate of jet production processes is too high to allow for a collection of all these events: the 
jet triggers are therefore prescaled, such that only a limited integrated luminosity is collected 
by them. The best dataset for a search of Z ^ bb decays is instead the one collected by 
a trigger requiring the presence of a clean muon candidate: the 6-quark jets produced in 
Z decays contain lots of low Pt muons, originated in B hadron decays, sequential charmed 
hadron decays, or other processes in smaller amounts. The starting point of the search 
was therefore the sample of 5.5 million events featuring a Pt > 7.5GeV/c central muon 
candidate, corresponding to an integrated luminosity of L = 103 ± 7 pb~^. 

2.1 The Datasets 

The initial dataset was cleaned up by requiring the muon candidate be also identified as 
a good muon by a standard offline filter. Jets were reconstructed by a fixed R = 0.7 cone 
clustering algorithm^, and two jets in each event were required to contain charged tracks 
forming a well identified secondary vertex (hereafter "5* FX ta(f^ , from the name of the silicon 
vertex detector). This selection reduced drastically the dataset, being satisfied by only 5,479 
events. The selection discarded the non-heavy flavour component of the QCD background, 
and increased the signal/noise ratio by about two orders of magnitude. 

We used the PYTHIA Monte Carlo| to generate 1.7M Z ^ bb decays that were subjected 
to a detector simulation, flltered by a trigger simulation, and passed through the same offline 
selection used for the real data. Using the predicted cross section for the searched process, 
the number of Z ^ bb decays in the double SVX tagged dataset is expected to amount to 
124 ± 14 events. 

As will be explained in the following, we have used real Z decays to electron-positron 
pairs collected during run 1 to better understand the behavior of initial state radiation in 
the Z production. These events were selected from high Et electron triggers by requiring 
the presence of a good quality central electron plus another electron candidate passing 
looser cuts. The dataset consists in more than six thousand events, and can be useful for 
kinematical studies, particularly when the variables one is interested in are difflcult to model 
with Monte Carlo programs. 

To study the contamination of other boson signals in our dataset we also generated 
500,000 W boson decays to cs pairs and 500,000 Z decays to c-quark pairs with PYTHIA 5.7: 
in fact, these processes may give a contribution to the muon dataset, due to the presence 
of charm quarks in the final state. These events were subjected to the same treatment 
described for the Z ^ bb data. Their contamination to the double SVX tags dataset was 
found to be totally negligible. 

F.Abe et alii, Phys. Rev. D 45 (1992), 1448. 
® H.Bengtsson and T.J.Sjostrand, Comput. Phys. Commun. 46 (1987), 43. 
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Sample 


Run 1 


Z ^bb 


Eff. 


Events in 






PYTHIA 




103 pb-^ 


Initial 




1,673,000 




(110± 12)10^ 


Trigger 


5,414,755 








2 SVX tags 


5,479 


1,867 


21.3% 


124 ± 14 


A<l>i2 > 3 


1,684 


1,368 


73.2% 


91 ± 11 


Sg^T < 10 


588 




(50.0%) 


(45 ±21) 



Table 1: Number of events selected by each of the cuts described in the text for the ex- 
perimental data and for the Z ^ bb MC data, efficiency of each cut for the Z signal, and 
expected Z events in the total dataset. The signal size and the efficiency of the last cut 
(TizEt < lOGeV) are estimated using Z — > e+e~ data. 

2.2 Kinematical Tools 

Even after the very restrictive selection of events with two SVX tags, the signal/noise ratio is 
too small to allow for a clear identification of the Z ^ bb signal in the dijet mass distribution: 
other tools are needed to increase the discrimination of the electroweak production of 6-quark 
jets from the strong interaction. 

From a theoretical standpoint, one expects the two production processes to be different in 
many ways. The two initial state partons of the Z production have to be a quark-antiquark 
pair, and the process is time-like. On the contrary, in the QCD creation of a bb pair both a 
quark-antiquark and a gluon-gluon pair can give rise to a time-like direct production process, 
and also space-like diagrams may contribute. Many of the QCD processes are expected to 
result in a pair of outgoing partons with a flatter pseudorapidity distribution than those 
from the searched Z decay; but the double SVX tagging and the requirement of a central 
muon candidate in one of the two jets result in pseudorapidity distributions that are already 
very constrained and well peaked at zero, due to the acceptance of the SVX and the muon 
chambers. 

If one examines the color structure of the diagrams one however notices a marked dif- 
ference between Z ^ bb and g — > bb. In the QCD processes there is a color connection 
between the initial and the final state absent in the Z production. Furthermore, while both 
the initial and the flnal state of the Z production are in a color singlet configuration, the 
opposite is true for QCD tree level diagrams. These considerations alone may lead us to 
believe that QCD gives rise to processes with a higher color radiation accompanying the two 
6-quark jets; furthermore, the pattern of this radiation is different. In QCD processes, in 
fact, the color connection present at LO between the two final state partons and the initial 
state should give rise to a enhanced radiation fiow in the planes containing each of the two 
leading jets and the z axis: color coherence prescribes the radiation from the incoming and 
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outgoing partons to interfere constructively in these regions, while the color singlet produced 
in the Z decay will emit soft radiation mainly between the two leading jets. 

Although no event-by-event discrimination appears possible by the use of variables that 
try to pinpoint the differences in the radiation pattern, two variables that deal inclusively 
with these features prove useful for our search. These are the azimuthal angle between the 
two 6-quark jets, A$i2, and the sum of transverse energies of all the calorimeter clusters in 
the event beyond the two jets, Tj^Et'- they both have some discriminating power between 
a high radiation process and the colorless production of a dijet system; but, while the 
first one is easy to simulate for Monte Carlo programs, being relatively independent from 
the modeling of the underlying event, the second is critically dependent on the detailed 
features of the initial state radiation mechanisms. For an homogeneous comparison one is 
therefore bound to use the experimental data to understand the behavior of the T.^Ej'- The 
distribution of the SVX data (to be considered a pure background sample, due to the very 
low S/N ratio) in the plane of the two kinematic variables can thus be compared to that 
of experimental Z decays to e~^e~ pairs^, as shown in fig.|l[ The difference in the radiation 
flow for these two samples is evident. 



CDF PRELIMINARY 




Inclusive Muon Data Z — > e"^ e' Data 



Figure 1: These plots prove the discriminating power between a QCD and a EWK process 
of the two kinematic variables we are selecting our data with. On the left is shown our SVX 
tagged data, on the right the Z — >■ e"'"e~ data. The two distributions are normalized to equal 
volume. 

We therefore select our data by placing tight cuts on these two variables: we require the 

^ These two variables have been shown to be rather insensitive to the final state radiation, that is present 
in the searched process but is absent in the leptonic final state of a Z decay. 
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two leading jets to be separated in azimuth by at least three radians, and the additional 
clusterized Et be less than 10 GeV. These cuts will be shown to maximize the expected 
signal significance in Section |3.1| , where all the ingredients necessary for the computation 
are presented. 



3 The Counting Experiment 

After the kinematic selection we expect the dataset to be still rich of background, with a 
S/N ratio as high as 1/6 for events with Mjj ~ 90 GeV/ c^. Under such circumstances, the 
mass distribution can be used to demonstrate the presence of a signal only if the background 
shape is very well under control. 

We use events where only one of the two leading jets carries a SVX tag as a background- 
enriched sample, and look for a signal as an excess of events in the double SVX data; to 
obtain a background estimate we rely on the observation that the probability of finding a 
secondary vertex in a jet is independent on the value of the kinematic variables we have 
selected our data with. We divide our data into four subsamples: events that fall in the 
kinematic region selected by our cuts (A<l>i2 > 3 radians, T^^Et < 10 GeV, in the following 
referred to as the Signal Zone"), and events that fail those requirements, from both the 
double tagged dataset — "(++)" in the following — and from events having only one tagged 
jet — hereafter "(+0)" events. We then define a tag probability as the ratio between (++) 
and (+0) events outside the Signal Zone, and extrapolate it inside, obtaining an absolute 
background prediction for the double tags in the Signal Zone: 

exp,in obs.in \ obs,out/ obs.outJ' 

The procedure just outlined is carried out for each bin of invariant mass of the dijet system: 
by doing that, we obtain a tag probability as a function of the dijet mass, from which 
we evaluate an absolute background prediction iVj^p j„(Mjj) that can be compared to the 
observed count in each bin of the mass distribution. An excess around 90 GeV/ (? will be 
bias-free evidence for the signal, if away from there the excesses are compatible with zero. 



3.1 The Choice of Cuts 

In order to extract the highest possible significance from the excess of events in the dataset 
we have studied the expected signal significance as a function of the cuts on the kinematic 
variables used for the selection. Our definition of the significance for the present purpose is 
an approximation: we define it as 5 = ISl signal/ {ptgr ® ^hgr ® \/^tot)- The Monte Carlo is 
used to estimate Nsignai as a function of the cuts, while Nfot is defined as the sum of Ngignai 
and the expected background, computed with the method described above. The statistical 
and systematic errors on the background prediction are added in quadrature to the Poisson 
fiuctuation of the total expected number of events Ntot, giving a value of S per each choice 
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of the cut on the variable studied. We are thus able to decide what are the best possible 
cuts on the two selection variables. 
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£,E|. cui definiiig ihe Signal Zone (GeV) A*,, cm defining Ihe Signal Zone (radians) Half Bin Size around (GeV/c"] 

Figure 2: Maximization of the expected signal significance as a function of the cut on 
the kinematic variables. Left: choice of cut on Tj-^Et; center: choice of the cut on A$; 
right: choice of the bin size. Four different possibilities for a systematic uncertainty in the 
background evaluation are considered. 

As can be seen in fig.0, the chosen cuts on A$ and Tj^Et do maximize the expected signal 
significance S, regardless of the systematic uncertainty attributed to the extrapolation of 
the tagging probability. 

The same machinery can be used to decide the optimal binning in the dijet mass distri- 
bution for the counting experiment. The dijet mass is of course a very discriminant variable: 
by defining appropriately the width of the bin around 90 GeV/ (? where we look for an excess 
of events over the background, we can again maximize the expected significance. Using the 
background prediction and the expected shape of the signal mass peak from Monte Carlo, 
we obtain significance curves that point to the best binning. The latter is chosen to be 
40 GeV I (? wide, as shown in fig.|[ 

3.2 Results of the Counting Experiment 

Having tuned cuts and binning to their most favorable values, we can perform the compu- 
tation of the background and study the results. These are detailed in table ^ 

The excess in the third bin is quite significant: using a conservative estimate of the 
systematic uncertainty in the extrapolation of the tag probability (a total of 4% is estimated 
from comparisons of observed and predicted events in signal-depleted samples) the prediction 
becomes Nexp = 248.49 ± 13.38 events. The probability of a fiuctuation of this number to 
the 318 observed events is 0.00061, equivalent to 3.23 standard deviations for a one-sided 
gaussian distribution. 
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Figure 3: Results of the counting experiment with a 10 GeV/c^ binning. The excess 
over background predictions is compared to the expected shape of a Z ^ bb mass signal 
(PYTHIA). The inset shows the mass distribution of double tags (points with error bars) 
and the expected background (full histogram) . 

The numbers quoted above are already compelling evidence for the presence oi Z ^ hb 
decays in our dataset. But we can use a smaller bin size to study the shape of the excess 
of events in the mass distribution: we then expect to see a gaussian peak0 at QOGeV^/c^, 
with a r.m.s. of 12.?)GeV / (? . The results of the counting experiment with this smaller bin 
size confirm the expectations: the excess fits very well to the expected signal shape (fig.^. 

We have still another method to verify the assumption that we are observing a Z signal. 
In fact, the counting method allows a study of the behavior of the excess as a function of 
the same kinematic variables used for the data selection, that have been shown to have a 
distinctive shape for the electroweak process. To do that, we select (++) and (+0) events in 
the Signal Zone that fall in the interval 70 < Mjj < 110 GeV/c'^: as table ^ shows, we have 
an excess of 69.5 ± 20.0 events there. We can then build the T^^Et and A$ distributions for 

A tuning of the jet energy corrections, appropriate for 6-quark jets where one of the two undergoes 
a semileptonic decay to a muon, has been performed to increase the evidence for the signal in the mass 
distribution. It aUowed us to increase the relative mass resolution aM/Mjj by 50% and to bring the 
reconstructed dijet mass to the nominal value in Monte Carlo events. 
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Mass Interval 


^JDSfcilVfciQ 




Il/XLt;SS 


Il/Xp. Zj 


0-30 GeV/c^ 
30 - 70 GeV/(? 
70 - llQGeV/c^ 
110- 150Ge1//c2 
150 - 190 GeV/c^ 
190 - 230 GeV/c^ 
230 - 270 GeV/c^ 
270 - 310Gey/c2 
310 - 350 GeV/(? 
350 - 390 GeV/c^ 


n 
u 


n nc; -u n no 
u.uo It u.uy 


n nc; -u i /in 
— u.uo It i.4U 




163 


149.30 ±6.13 


13.70 ± 14.16 


1.7 ±0.8 


318 


248.49 ± 8.95 


69.51 ± 19.95 


41.4 ± 18.8 


66 


65.51 ±4.68 


0.49 ±9.38 


2.4 ± 1.1 


29 


19.70 ±2.59 


9.30 ±5.97 




7 


7.06 ± 1.73 


-0.06 ±3.16 




3 


2.24 ±0.82 


0.78 ±2.38 




1 


1.60 ±0.95 


-0.60 ±2.00 







0.55 ±0.98 


-0.55 ± 1.71 




1 


0.00 ±0.00 


1.00 ± 1.76 





Table 2: Results of the counting experiment. 
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Figure 4: The excess observed in the 70^110 GeV bin is distributed according to expectations 
in the T^^Et variable (left) and in the A$ variable (right). 



the (+±) events and compare them to the corresponding (±0) distributions scaled down by 
the tag probability P90 = iV„"^+o„i(70 ^ 110)/iV„"^°o„j(70 ^ 110). Indeed, the double SVX tags 
in excess are distributed in the two kinematic variables as expected for Z events, as can be 
seen in fig.^. 
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4 Unbinned Likelihood Fits to the Mass Distribution 



Having established the presence of a. Z bb signal in our dataset with a counting method, 
we can perform a two-component fit to the dijet mass distribution, and extract additional 
information from its shape. 

We used a three-step procedure to fit the data. First of all, we obtained a background 
shape from a unbinned likelihood fit to a signal-depleted sample consisting in (-1-0) events 
falling outside of the Signal Zone. This sample is expected to be totally dominated by the 
QCD background, and can be successfully fit to the following simple functional form: 

PbgriMn) = i-e-^-/^®iV(/i,a) 

where "®" is the convolution operator, while N{fi, a) is the Normal distribution. The fit is 
shown in the lower plot of fig.^. 

The above form cannot be directly used as a background shape for double SVX events, 
since the probability of tagging a second jet is correlated to the invariant mass. To account 
for that, we obtained a tag probability curve P{Mjj) by performing a fit to the ratio of 
double and single SVX tags rejected by the kinematic cuts. 

The knowledge of the parameters of the two above fits allows us to write a two-component 
unbinned likelihood for the dijet mass distribution of the sample, as follows: 



N++ 



n 



sig 



■Ps.,(M.)+n+;-P4,(M,) 



"^sig ~l~ '^bgr 



where: A^++ is the observed number of events in the {++) sample; and are 

respectively the number of signal and background events in the sample; Plp.g{M) is 

the normalized background mass shape of the events in the (-I--I-) sample, obtained from the 
background fit to the single SVX data and the tag probability curve; and, finally, Psig{M) 
is a normalized gaussian describing the signal. 

By maximizing log C with respect to the number of signal and background events, and to 
the mass and width of the signal shape, but keeping frozen the background shape parameters 
and the tag probability curve parameters obtained previously, we obtained = 91 ± 30 
events, Mz = 90.0 ± 2.4 GeV^/c^, and az = 9.4 ± 3.5 GeV^/c^; the fit results are shown in the 
upper plot of fig.^. The mass and width of the gaussian peak are in perfect agreement with 
Monte Carlo expectations (M = 90.0 GeV/c^, aM = 12.3 GeV/c^), while its normalization 



is in agreement to the excess obtained in Section but is slightly larger than the Monte 
Carlo predictions {Nz = 45.5 ± 20.7 events). Although partially correlated to the results 
of the counting experiment, these numbers give additional qualitative informations on the 
behavior of the signal we have isolated. 

Various studies of the possible systematic uncertainties of the fitting procedure have 
been performed. The systematic uncertainty due to the negligence of a signal contamina- 
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Figure 5: Top: results of the unhinned likelihood fit to the dijet mass distribution in the 
(++) sample; bottom: results of the unbinned likelihood fit to the dijet mass distribution of 
(+0) events rejected by the kinematic cuts. 



tion in the background sample was estimated to be ±9 events; the systematics due to the 
parametrization of the tag probabihty amount to ±16 events; and, finally, the systematic 
uncertainty due to the modeling of the background shape was showed to amount to ±5 
events. Therefore the final result of the unbinned likelihood fit to the mass distribution is 
the following: A^'^ = 91 ± 30(stat) ± l%syst) events. 



5 Conclusions 

We have searched 5.5 million inclusive muon events collected by CDF during run 1 for a 
Z — > 66 signal. The very low signal/noise ratio at trigger level (less than 10~^) implies 
that a really strict selection is required in order to isolate a signal. We have designed an 
optimized selection for that purpose by making use of double SECVTX tagging plus some 
kinematic criteria suited to the discrimination of an electroweak production from the QCD 
background. By these means we select 588 events, 318 of whose have a reconstructed dijet 
invariant mass between 70 and 110 GeV/(?, where the Z decay is expected to yield 41.4±18.8 
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events. By comparing the observed events with the expected background we find evidence 
of the signal, quantifiable in an excess of 69.5 ± 20.0 events having a suggestive shape in 
the mass distribution. The excess corresponds to a 3.23cr fiuctuation of the background, 
in the hypothesis of no signal. Finally, we have used a unbinned likelihood fit to the dijet 
mass distribution to obtain additional evidence of the presence of 91 ± 30 (stat) ±19 (syst) 
Z ^ bb decays in the signal sample. 
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